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Abstract 

Electric field is a thermodynamic force that can drive collisional inter-ion-species transport in 
a multicomponent plasma. In an inertial confinement fusion (ICF) capsule, such transport causes 
fuel ion separation even with a target initially prepared to have equal number densities for the two 
fuel ion species. Unlike the baro-diffusion driven by ion pressure gradient and the thermo-diffusion 
driven by ion and electron temperature gradients, electro-diffusion has a critical dependence on 
the charge-to-mass ratio of the ion species. Specifically, it is shown here that electro-diffusion 
vanishes if the ion species have the same charge-to-mass ratio. An explicit expression for the 
electro-diffusion ratio is obtained and used to investigate the relative importance of electro- and 
baro-diffusion mechanisms. In particular, it is found that electro-diffusion reinforces baro-diffusion 
in the deuterium and tritium mix, but tends to cancel it in the deuterium and helium-3 mix. 
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I. INTRODUCTION 



In inertial confinement fusion (ICF) experiments, where the fuel assembly is a binary 
mixture of deuterium and tritium or deuterium and helium-3, the fusion power production 
is proportional to the product of number densities of the two species. Optimal fusion yield 
at a given target temperature in local thermodynamic equilibrium requires not only the fuel 
assembly be equi-molar for the two ion species, but also of equal number densities everywhere 
in the assembly [1]. In terms of the mass densities, this suggests that the mass density of 
the light ions 

pi = mini 

should be mi/m^ times of the heavy ion mass density 

p h = m h n h 

for n\ = rih- Here the subscript "1" denotes light ion species and "h" denotes heavy ion 
species. Defining the mass concentration of the light ions as 

c = — 

~~ P 

with the mixture ion mass density 

P = Pi + Ph, 

one finds that the optimal arrangement of equal ion number densities implies a spatially 
uniform c = mi /(mi + rrih)- This condition can be accurately satisfied in the initial target 
preparation. The dynamical process of implosion, however, can introduce light and heavy 
ion separation which degrades the fusion power production [2-4]. 

The collisional inter-ion-species transport or concentration diffusion is driven by the con- 
centration gradient Vc as well as other thermodynamic forces such as the ion pressure 
gradient, electron and ion temperature gradients, and electric field. In the case of a plasma 
with two species of ions, we show that the diffusive ion mass flux i takes the general form of 

i = - p D(y C + k p W log ^ + + kf V logT; + 4 e) V log T e ) . 

This flux governs c evolution through 

dc 

p— +pu- Vc + V-i = 0, (1) 
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where u is the plasma fluid velocity. 

For the ICF fuel assembly, even when the initial condition has Vc = by design, the 
theromodynamic cross terms can drive significant diffusive flux i through the baro-diffusion 
(k p 7^ 0), electro-diffusion (fc^ ^ 0), and thermo- diffusion (fcy 7^ 0)- The above equation 
attributes baro-diffusion to the total ion pressure gradient Vpj, which is the sum of the 
ion species pressure gradients. This is different from the case of the neutral gas mixture, 
where baro-diffusion is considered to be due to the total mixture pressure. Also, because 
of the large difference between the electron and ion masses their temperatures can vary. 
In the general case, thermo-diffusion is driven by gradients of both the electron and ion 
temperatures. 

Baro-diffusion is fundamentally the result of the mass depedence of the species thermal 
speed, and we will show it to be indepedent of electric field and effective gravity due to accel- 
eration and decceleration in implosions of the ICF target, just like its counterpart in neutral 
gas mixtures [5, 6]. This finding should be contrasted with a previous result [2, 3] which 
suggests gravity and electric field dependence of k p for binary plasma mixture. The role of 
electric field, or electro-diffusion, is a feature intrinsic to a plasma. It is fundamentally the 
result of the different acceleration experienced by the ions of different charge-to-mass ratio 
in an electric field. In a low temperature plasma with significant neutral gas background, 
this effect is known as ions having different mobilities [7]. We show that electro-diffusion 
vanishes (fc^ = 0) if the charge-to-mass ratio is identical for the two otherwise distinct 
ion species. The charge-to- mass ratio dependence of &e allows drastically different electro- 
diffusion behavior for different binary plasma mixtures. For the DT mix, we find that 

fee = k p , 

so the electro- and baro-diffusion reinforce each other in a plasma shock. In contrast, for 
the D 3 He mix, one has 

Ke — ~kp, 

which implies that electro-diffusion tends to cancel baro-diffusion in a plasma shock. 

The thermo-diffusion is fundamentally the result of the thermal force in the collisional 
drag between different ion species and between electrons and ions. The ion-ion thermal 
force produces a non- vanishing kj), while the ion-electron thermal force leads to a finite 
fcjf ' . Because of these dependences, thermo-diffusion coefficients require a kinetic calculation 
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of the said transport coefficients through the distribution function perturbed from a local 
Maxwellian. Hence, are not thermodynamic quantities like k p and /eg. It must be noted 
that with comparable ion masses, the thermal force evaluation is a much more involved 
exercise, than the electron-ion one carried out by Braginskii [8]. The critical information 
is nevertheless implicitly contained in standard transport calculations for multi-component 
plasmas such as that by Hirshman and Sigmar [9]. Explicit evaluation of the thermo-diffusion 
coefficients for the DT and D 3 He mixtures will be carried out in a future work. 

The inter-ion-species diffusion, which modifies the relative number density of the two fu- 
sion reactants, has been an issue of interest in inertial confinement fusion and dense plasma 
research. Among those we are familiar with, C. H. Chang, B. Albright, and W. Daughton 
from Los Alamos National Laboratory have investigated the models for i with varying de- 
grees of approximation for evaluating the baro- and thermo-diffusion in the past decade, 
albeit in unpublished reports. More recently, it has been attracting special attention for the 
pioneering analysis of Amendt et al [2, 3], which shows baro-diffusion may be responsible 
for the discrepancy between the neutron yield measured during ICF implosions and that 
predicted by simulations. This can largely be attributed to the realization that the strong 
pressure and temperature gradients, along with the strong electric field [10], are induced 
by shock waves inevitably present in the imploded capsule. Inter-diffusion between the two 
ion species must therefore take place; the resulting separation of the fuel constituents in 
the hot spot can significantly degrade the fusion yield [2, 3]. An experimental evidence for 
the fuel stratification in ICF implosions has been reported by Casey et al based on their 
measurements at the OMEGA laser facility [11]. 

The purpose of this paper is to clarify the underlying diffusion mechanisms, especially 
the role of electro-diffusion in relation to baro-diffusion. The derivation of i based on multi- 
component collisional fluid models also provides a framework for incorporating this impor- 
tant physics in ICF modeling. As highlighted earlier for a binary plasma mix, there are a 
number of subtleties absent in neutral gas mixture, and not addressed in preceding studies. 
Since the inter-species diffusion is a general topic for multi-component plasmas, it is ex- 
pected that the results described here would be of general interest, for example, to tokamak 
edge plasma modeling, plasma processing, and stellar structures. 

The rest of the paper is organized as follows. In the next section, the thermodynamic 
framework for evaluating electro-diffusion is outlined. In section III, the momentum equation 
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for an ion species in the center-of-mass frame is obtained and simplified by imposing an 
ordering relevant to a collisional plasma shock. Then, in section IV, this equation is utilized 
to write the diffusive mass flux of an ion species in the general form of section II, thereby 
explicitly evaluating the electro-diffusion ratio. Finally, in section V, implications of electro- 
diffusion in the ICF context are discussed. 



II. THERMODYNAMIC EXPRESSION FOR THE DIFFUSIVE FLUX 

To define the diffusive mass flux of ion component a, the center-of-mass velocity u is first 
introduced by 

pU = J2pa U a^ ( 2 ) 

a 

where p a is the partial density of component a and p = J2 a Pa is the mixture total ion 
density. Also, u a denotes the net flow velocity of the component and the sum on the right 
side of Eq. (2) is over all the components present in the mix. The diffusive mass flux is then 
given by 

i Q = p Q (u a -u). (3) 

If the system is close to local thermodynamic equilibrium, a linear relation between the 
thermodynamic forces and the resulting fluxes can be assumed. For a neutral gas mixture, 
the total diffusive mass flux of component a can be written as [5, 6] 

i Q = -pD ( Vc Q + k p V log p + k T V log T) , (4) 

where p and T denote the mixture total pressure and temperature, respectively, and c a = 
p a / p denotes concentration of the component a. Parameter D is called diffusion coefficient; 
it governs the diffusive flux when only the concentration gradient is present. In view of 
Eq. (4), baro- and thermo-diffusion coefficients are then equal to k p D and /cy-D, respectively. 
Dimensionless parameters k p and kr are usually referred to as baro- and thermo-diffusion 
ratios, respectively. 

Interestingly, k p is a thermodynamic quantity, i.e. it can be evaluated given local values of 
thermodynamic variables and does not depend on the details of collisions [5] . In particular, 
for a binary mix it can be found [5, 6] 

k p = c(l-c)(m h -m l )( K ^ + ^—^), (5) 
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where mi and rrih are the atomic masses of the light and heavy fractions, respectively. 
Also, c = q stands for the concentration of the light fraction; the subscript "1" is dropped to 
simplify notation, as in a binary mix the concentration of the heavy fraction can be recovered 
through 0^ = 1 — q. In contrast to k p , k? is an intrinsically kinetic quantity and is subject 
to change depending on the collisional model. 

To account for the effect of the electric field, the corresponding force needs to be added 
on the right side of Eq. (4) to rewrite it as 

i a = - P D(Vc a + fcpVlogp + fcrVlogT + ^V$), (6) 

where $ is the electrostatic potential and, by analogy with k p and fcr, the electro-diffusion 
ratio &e is introduced. In what follows, we focus on the case of a binary plasma mix, i.e. 
plasma consisting of two ion species and electrons. To illuminate the new features brought by 
the electric field, as compared to the case of a binary mix of neutral gases, we apply Eq. (6) 
to the system consisting of the two ion species. Within such an approach, the electron 
species is viewed as an external factor that affects the system of interest through the electric 
field and collisions. In other words, the electrons do not make contribution in the definition 
of p, u, and p, that is consistent with our objective to understand the relative motion of the 
two ion species. It is worth noticing that because of small inertia the electron contribution 
to the overall plasma density and flow is negligible. Thus, for all practical purposes p and 
u can still be referred to as the plasma density and flow, respectively. On the contrary, the 
electron and ion pressures are generally comparable and employing the total ion pressure, 
rather than the overall plasma pressure, in place of p in Eq. (6) is crucial. 

Assuming the diffusive flux of the form (6) it is possible to evaluate \ze by generalizing 
formal thermodynamic methods used in Ref. [5] to evaluate k p [4]. Instead, here we start 
from the first-principle based momentum conservation equations for individual species to 
automatically recover this form. Importantly, in addition to readily providing k p and ks, 
this technique gives formulae for D, kj) and k^ in terms of standard transport coefficients. 
In so doing, it lays the framework for evaluating the overall effect of the ion concentration 
diffusion that is inherently not possible within the thermodynamic approach. 
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III. MOMENTUM CONSERVATION FOR ION SPECIES 



We start by writing the momentum equations for the two ion species 

where the subscript a can be "1" and "h" to denote the light and heavy ion species, respec- 
tively In Eq. (7) n a , p a and Z a stand for the species' number density, partial pressure and 
charge number, respectively. Also, E stands for the electric field. Acceleration due to an 
external force of a non-electric origin, such as the gravitational force, is denoted by F Q , while 
Rq,^ is the force density due to collisional momentum exchange with the species (3 and the 
sum on the right side of Eq. (7) is over all plasma species, including electrons. The pressure 
tensor P a is defined in the frame co-moving with the species net flow by 

P a = m a / (v - u a )(v - u a )f a d v, 

where f a is the species distribution function. Finally, 

d a /dt = d/dt + u Q ■ V. 

Next, we obtain the momentum equation for the center of mass velocity defined by Eq. (2). 
To do so, it is convenient to rewrite Eq. (7) in the conservative form: 

" + V • (p a u a u a ) + V • P a - n a Z a eE - p a F a = ^ R Q/3 , (8) 
at p 

where species continuity equation 

daPa/dt + p Q V ■ u Q = 
is used. Introducing the species' velocity in the center-of-mass frame 

w Q = u Q - u, (9) 

noticing that 

V • (p a u a u a ) = V • (pqUu) + V • (p a w a u) + V • (p a uw a ) + V • (p a w Q w Q ) 
and summing Eq. (8) over all the ion species we find 

+ V • (pan) + V • Pi - XJ (n a Z a eE + p a F a ) = £ (10) 

a=ii a=ij 



where we use that Yl a =i p Q w Q = 0, Pi = J2 a (Pa + |O a w a w tt ) is the total ion pressure tensor 
in the center-of-mass frame and subscript "a = if denotes summation over the ion species 
only. Equation (10) is then easy to transform to a more familiar form 

P ^dJ + VPi + V ' " E (n a Z a eE + p a F a + R Qe ) = 0, (11) 



where d/dt = d/dt + u-V and pi and flj are the total ion pressure and viscous stress tensor, 
respectively 

Finally, we rewrite equations for the individual species flows in the center-of-mass frame 
to obtain 

dpa ™ a + Pa w a ■ Vu + p a w Q V ■ u + V ■ P c a m - n a Z a eE - p a F a + p a ^ = J2 R «/3' ( 12 ) 
dt dt p 

where 



P a + p Q w a w a = m a y(v - u)(v - u)f a d 3 v 



is the species' pressure tensor in the center-of-mass frame. By splitting a scalar pressure out 
of P^,' m , Eq. (12) is then rewritten further to find 

^^+ Pa w a -Vu+p a w a V-u+Vp c a m +V-U c a m -n a Z a eE-p a F a +p a ^ = ^R^, (13) 
ere dt p 

where and n^* m are the species' partial pressure and viscous stress tensor, respectively, 
evaluated in the center-of-mass frame. 

The calculation presented here is motivated by the problem of the ion species diffusion 
within an ICF relevant shock wave front. In general, the shock front width can be as small 
as the mean free path, making a local treatment, as well as the framework of the previous 
section, invalid. However, for a moderately strong shock, the front width can be assumed 
much greater than the mean free path A, i.e. 

A/A « 1, (14) 

where A is the characteristic spatial scale of the plasma (e.g. the shock width). The mean 
free paths can substantially differ for pre- and post-shocked plasmas; for defmitiveness, we 
refer A to the post-shock mean free path. Ordering (14) is usually satisfied for Mach numbers 
< 2 and ensures that the plasma remains mostly collisional throughout the shock front [12]. 
The characteristic temporal scale r can then be estimated from r _1 ~ v s h/A, where v s h 
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is the shock speed. Assuming that ion masses are comparable across different species and 
introducing the characteristic ion thermal speed ikh-i ~ v s h we find 

Vl r~ (A/A)- 1 » 1, (15) 

where z/j ~ vth-i/^ is the characteristic ion collision frequency. 

Estimating the friction between the ion species by nihUiiiiiwi — Wh), where /i^ is the 
reduced mass for the light and heavy ions, it is straightforward to show that ordering (15) 
makes 

A « 1, 



vth-i A 

thereby ensuring that the system is close to a local equilibrium. The terms on the left side of 
Eq. (13) that contain both the spatial gradient and w a are quadratic in the small parameter 
A/A and can be dropped. For the same reason P c a m ~ P a and the superscript "cm." 
appearing next to the partial pressure and viscous tensor can be omitted. Finally, due to 
the same estimate for the friction, the dw a /dt term can be dropped as well and Eq. (13) 
reduces to 

<-> du x , 

Vp a + V ■ Il a - n a Z a eE - p a F a + p a — = ^ Ra/3- (16) 

dt p 

Equation (16) is valid for a plasma with an arbitrary number of species as long as order- 
ing (15) is obeyed. In the next section, we apply it to evaluate the electro-diffusion coefficient 
in a plasma with two ion species. 



IV. EVALUATING ELECTRO-DIFFUSION COEFFICIENT 

The viscous term appearing on the left side of Eq. (16) is governed by the second order 
derivatives of macroscopic parameters. Its contribution is therefore not retained in Eq. (6), 
which is obtained by assuming linear relation between the thermodynamic forces and the 
resulting flux. In principle, this contribution may be substantial and effectively modify the 
baro-diffusion ratio [5, 6]. However, as the main goal of the present study is to elucidate the 
role of the electric field on the ion diffusion, in what follows we drop V • H a on the left side 
of Eq. (16). Then, employing Eq. (11) with the V ■ IL, term also dropped to evaluate du/dt 
in Eq. (16) we find 
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(Vp Q - ^V Pi ) - (Z a n a - ™ ]T Z p n p )eE - (p a F a - ^ £ ppFp). (17) 

P P fS =ij P =ij 

Notice, that if we had included electrons into the system, the sum over (5 in the second term 
on the right side of Eq. (17) would vanish due to quasi- neutrality. 

Equation (17) gives the light ion species diffusion velocity w/ through the Ry, dependence 
on the net velocity difference between the ion species since 

- w h = (l + — )wi = - — — , (18) 
v ph J 1-c 

where c is the concentration of the light ion species. In a multi-component plasma [9], the 
collisional drag between species a and j3 

a 

VT a + c%npVT P l (19) 

where, in general, coefficients A a p, c^l and c^ 2 j are complicated functions of the masses, 
densities and charge numbers of all the species and H aa = for arbitrary T a implies = 
c ai = 0- Also, T a is the temperature of species a, fi a p = m a mp/(m a + mp) is the reduced 
mass and v a p stands for the frequency of collisions between species a and (3. Conventionally, 
the terms on the right side of Eq. (19) proportional to the velocity difference and temperature 
gradients are referred to as the frictional and thermal forces, respectively. Summing Eq. (19) 
over the ion species we find 

R-q/3 = - [A a pp a pn a v a p(v? a - wp) + B a pnpVTp], (20) 

where B aa = Y.p=i 3 c^p and B a p = c^j for a ^ (3. 

When the elementary masses of species a and /3 are comparable, the thermal force acting 
between them depends on both VT Q and VTp, i.e. c£i ~ c^ 2 j on the right side of Eq. (19). 
In contrast, the thermal force acting between the electron and any of the ion species is 
dominated by the electron temperature gradient, because the thermal speed of electrons is 
much greater than that of ions, i.e., 

R ae = -{A ae p ae n a v ae (w a - w e ) + B ae n e VT e ] (21) 

where we set B ae = c^} to unify notation with Eq. (20), and 

w e = u e - u 
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is the electron flow velocity in the ion center-of-mass flow frame. Unlike ~ (X/A)vth-i 
from our collisional ordering, there is no such constraint on the electron flow w e due to the 
much larger electron thermal velocity. This is consistent with the well-known result that a 
collisional plasma can carry a substantial current in the electron channel despite that the ion 
current is negligibly small in the short mean-free-path limit. To estimate the ion-electron 
frictional force on the right side of Eq. (21) we introduce the plasma current 

J = — en e u e + ^2 en a Z a u a = — en e w e + ^ en a Z a w a , (22) 

where the quasi-neutrality condition along with Eq. (9) is used to obtain the right side of 
the equation. Then, using Eq. (18), the friction between the light ions and electrons can be 
rewritten as 

R>f e = -A le fii e mu le (wi - w e ) = -Ai e fii e mu le \—- r7T~~^w / w i + ~ J 1- ( 23 ) 

l cZi/mi + (1 - c)Z h /m h en e 

The term on the right side involving is smaller than the friction between the ion species 



by a factor of ym e /m^h. In an ambipolar plasma J = and the frictional force can be 
neglected on the right side of Eq. (21). Moreover, even for a plasma carrying significant 
current through the electrons due to the — en e w e term in Eq. (22), the ion-electron friction 
force is much less than its ion-ion counterpart as long as 



3 (24) 



en e v th _i Ay m e 

Condition (24) is the most restrictive in the case of a weak shock, where the shock front 



width can be many times of the ion-ion mean free path making (A/A)y m^h/m e of order 
unity [12]. The constraint on the plasma current becomes J <C en e v t h-i- Hence, in the 
absence of large currents on the order of en e v t h-i or greater, R^ e can be ignored for an ICF 
relevant shock wave. Consequently, the ion-electron collisional drag is dominated by the 
thermal force, 

R« e « -B ae n e VT e . (25) 

Applying general expressions (20) and (25) to our case and setting a = I the left side of 
Eq. (17) is now evaluated to find 

^ R«/3 + (Rie — — X! ^/3e) = 

0=ij p p=i 3 

-Ai h fii h mui h (wi - w h ) - BipnpVTp - (B ie - c ^ B Pe )n e VT e . (26) 

fi=l,h 0=l,h 
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While the coefficients Bp e are relatively easy to recover by generalizing the corresponding 
Braginskii's result for a simple plasma [8], evaluating Aih, B\h and Bu is quite complicated 
even in the case of only two different ion species with comparable masses, charge numbers 
and concentrations. Fortunately, it will be found unnecessary for the purpose of this paper, 
so we proceed leaving coefficients of Eq. (26) unspecified. 

To complete the calculation, Eq. (17) needs to be rewritten in the canonical form (6). 
The terms on the right side of Eq. (17) then have to be expressed in terms of the total ion 
pressure pi and the light species concentration c. We now proceed by doing so in the first 
term on the right side of Eq. (17) that is responsible for baro-diffusion. 

First, we observe that energy exchange between the ion species with comparable masses 
takes place over the same time scale as thermal equilibration within any of the two species. 
Hence, under ordering (15), TJ « Th and the overall ion temperature Tj can be introduced. 
Next, we notice that Pi = (nj + n h )Ti = (pi/mi + ph/m h )Ti, where m ; and are the light 
and heavy ion masses, respectively, to obtain 

Pi = -Tj n — Pi, 27 

cm h + (1 — c)mi 

and 

„ cmh _ Pirrihmi _ . . 

V ^ = Z7T S + r V n \ 12 Vc - 28 

cm h + (1 - c)mi [cm h + (1 - c)mi\ A 

The first term on the right side of Eq. (17) is then evaluated to find 

(Vpi - ^V Pl ) = — [Vc + c(l - c)(m h - mM— + — -)viogpi]. (29) 

Expression inside the square brackets of Eq. (29) is normalized, i.e. the coefficient in 
front of the Vc term is equal to unity. In view of Eq. (6) it means that the coefficient in 
front of the Vlogp^ term is equal to the baro-diffusion ratio k p . Importantly, this ratio, 
obtained here from ion fluid equations, matches the result (5), found in Refs. [5, 6] for a 
binary mix of ideal gases. Of course, this is just a reflection of the aforementioned fact 
t licit Jvp is a thermodynamic quantity and does not depend on the details of the collisional 
exchange between the species. It should be noted that recovering the same k p as in Refs. [5, 6] 
manifests the key difference between our approach and that of Refs. [2, 3], where k p is found 
to be dependent upon the electric field. 

With the technique presented in the preceding paragraphs, k^ can be straightforwardly 
calculated in the same way as k p . Before doing so, we take a brief detour and apply this 
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technique to clarify the role of gravity. The effective gravity appears in ICF relevant prob- 
lems when acceleration of the capsule during implosion needs to be accounted for. Upon 
switching to the frame co-moving with the capsule the inertial force enters the momentum 
equation that is formally equivalent to placing the system into external field with an effective 
gravitational acceleration g. 

In Refs. [2, 3], gravity is found to modify the expression for k p , so does the electric field. 
Within the framework of the present study, gravity can be included by setting the external 
force F a equal to g for both a = I and a = h. The third term on the right side of Eq. (17) 
is then found to vanish; that is, the gravitational force does not drive a diffusive flux. This 
result obtained with a rather formal method has a trivial physical explanation. Namely, 
gravity gives the same acceleration to all ions regardless of their mass and charge number 
and therefore introducing it into otherwise unchanged system does not directly contribute 
to the species separation. Of course, gravity can still affect ion concentrations indirectly. 
For example, it can do so by modifying the electron pressure balance. The electric field then 
has to adjust, thereby modifying the ion flux through its electro-diffusive component. 

Now we obtain the electro-diffusion ratio by writing the total diffusive flux of the light 
ion species in the canonical form. The second term on the right side of Eq. (17) is first 
evaluated to find 



(Z ini - ^ £ Zf,n p )eE = pc(l - c)(^- - — )eE. (30) 



Next, the terms on the right side of Eq. (17) are collected with the help of Eqs. (29) and (30) 
and Eq. (26) is employed along with Eq. (18) to find 

L = p lWl = -pD ( Vc + k p V \og Pi + ^ V$ + fc? V log Ti + k^V log T e ) , (31) 

where, as recovered by Eq. (29), k p is still given by Eq. (5) and 

D= - pT * x C(1 " C \ , (32) 
A lh p lh niui h cm h + (1 - c)mi 

k E = mim h c(l - c)(— H -)(— (33) 

v m; rrih ' x rni rrih' 



# = mi m h (— + —) rr« + ^ , (34) 

v m; rrih m i m h 

, CpI / c 1 — c\ ycZi (1 — c)Zhi r , , „ „ ,T f „ 
4 e) = rmmJ— + ) — l - + ± L-± 1 - c)B le - cB he ]^ } 35 

^m.i 177.1. ' L 177., 177.1. J / ■ 



/ c 1 — c 



cB u , (l-c)B lh 



-mi m h 



mi m h 



where quasi- neutrality condition was used to write Eq. (35). 
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Expression (31) does not have the exact form of Eq. (6) since the VlogT e term appears 
on the right side. This is because the only external force accounted for by Eq. (6) is the 
electric field, whereas for the system of the two ion species considered here the thermal force 
exerted by electrons is also external. Moreover, equation (6) is only valid when at any given 
point different components of the system are nearly equilibrated; in particular, this means 
that temperatures of all the components must be equal. For the system including ions only, 
this condition is satisfied due to our ordering (15). However, this ordering does allow T e to 
be different from T», as the energy exchange between the electron and any of the ion species 
takes longer than that between the two ion species by a factor of \jmi^Jm e . Hence, even for 
the plasma as a whole, for which the ion-electron thermal force is internal, T e and Tj have 
to be set equal for thermodynamically obtained Eq. (6) to be recovered. It is interesting to 
note that Tj ^ T e is normally expected in an ICF capsule, especially at the hot spot where 
fusion occurs. 

Equations (32) and (33) give the electro-diffusion coefficient feD, thereby fulfilling the 
goal of this paper. Notice that goes to zero if the charge-to-mass ratios are equal for the 
two ion species. This result rigorously obtained here from ion fluid equations has a simple 
physical explanation. Indeed, when Zi/rrii = Zh/rrih the electric field does not distinguish 
between the light and heavy ions and therefore does not contribute to the relative motion 
of the species. 

Unlike expressions (34) - (35) for thermo-diffusion ratios, which involve transport coeffi- 
cients B a/ 3, Eq. (33) provides an explicit result for fog without invoking a kinetic calculation. 
In other words, as its baro-diffusion counterpart, the electro-diffusion ratio is a thermody- 
namic quantity. Interestingly, it can then be evaluated in the same relatively simple way as 
suggested in Ref. [6] for evaluating the baro-diffusion ratio. We outline this calculation in 
the appendix A. 

V. DISCUSSION 

Relations (32) and (33) do not provide an explicit result for the electro-diffusion coefficient 
because of the transport coefficient A\h entering the formula for D. However, as the approach 
presented does provide an explicit result for fog, a substantial insight into the role of electro- 
diffusion can still be gained. To compare baro- and electro-diffusion caused perturbations 
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of the species concentrations we employ Eqs. (5) and (33) to evaluate the ratio 



k E _ Zi/mi - Z h /m h 
k p I /mi - l/m h 



(36) 



which depends on the properties of the ions only. In the special case of the two isotopes of 
one element, i.e. Z\ = = Z, Eq. (36) gives k E jk v = Z. In particular, for the practically 
important DT mix, the two coefficients turn out to be equal, k E /k p = 1. For the D 3 He mix, 
commonly used to study sub-ignited implosions, k E /k p = —1. In a plasma shock wave, the 
electric field is directed towards the unshocked region to prevent electrons' running ahead 
of ions and maintain quasi-neutrality. It can therefore be observed that in the case of the 
DT mix baro- and electro-diffusions act together, whereas in the case of the D 3 He mix the 
two tend to cancel each other. 

Of course, relation between the baro- and electro-diffusion ratios alone is not sufficient for 
relating the corresponding fluxes. The total ion pressure gradient and the electric field also 
need to be compared. While carrying out this comparison in a general case is hardly possible, 
it is reasonable to assume Vlog^j ~ eV$/Tj. Moreover, in a shock wave, the electric field is 
rather governed by the electron pressure gradient. The pressure of electrons is often greater 
than that of ions and therefore it is likely that the electro-diffusive flux may be noticeably 
larger than the baro-diffusive flux. This becomes particularly intriguing for the D 3 He mix, 
in which electro-diffusion counteracts baro-diffusion. As a result, 3 He concentration may be 
increased over its unperturbed value, contrasting the neutral theory based expectation that 
it is the lighter fraction whose concentration is enhanced in the shock front [13]. 

In terms of numerical modeling of the diffusive separation of the fuel ions in ICF capsules, 
the most direct approach would be to solve the mult i- component plasma equations in its 
individual species form. The electric field is then explicitly evolved. Alternatively, the ion 
fluid equations can be solved in the center of mass frame, i.e. p, u,pi, with the ion species 
concentration c followed by Eq. (1) and Eq. (31). With this approach, the electric field can 
be either independently evolved using Maxwell's equations, or inferred from the equation of 
motion for the electrons in the quasineutral regime. In this latter case, the electron inertia 
and electron viscosity are ignored, so 



The above equation implies that, at a minimum, the fluid equations should evolve the 




(37) 
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electron temperature separately from the ions', which fortunately is frequently done in ICF 
codes. 

Finally, we comment on whether or not electro-diffusion, described here by considering 
the ion species separately, can be attributed to baro-diffusion in the plasma as a whole. As 
previously mentioned, the plasma mass flux is essentially equal to the ion mass flux, because 
the electron inertia is negligible. The question to be answered is therefore whether or not 
the right side of Eq. (31) can be represented in terms of the total plasma pressure gradient, 
rather than in terms of the partial ion pressure gradient and the electric field. To investigate 
the issue, we insert Eq. (37) into Eq. (31) to obtain 

i, = - P DWc+ kp V Pi + \v Pe + *f Vlog7- + A;~ T (e) VlogT e j, (38) 
{ni + n h )Ti n e Ti 



where 



k T {e) = - (B le + B he ){T e /Ti)k E . (39) 



Eliminating both the electron and ion partial pressure gradients in Eq. (38) by substituting 
the total pressure gradient is only possible if 

^ = (40) 

for any values of n\ and nh- 

Combining Eq. (36) and the quasi-neutrality condition n e = Zini + Z^nh, one finds that 
Eq. (40) may be identically satisfied only for Z\ = Z^ = Z. The right side of Eq. (40) is then 
equal to Z and indeed matches the left side of Eq. (40) according to Eq. (36). Employing 
n e/(ni + rih) = k E /k p = Z in Eq. (38) we find 

= - p D [Vc + k p V logp + V log Ti + k T {e) V log T e ] , (41) 

where p = pi + p e is the total plasma pressure and 

kp = (1 + ZT e /Ti)k p . (42) 

Equation (36) predicts a larger baro-diffusion coefficient, as compared to the case of a neutral 
binary mix. In the limiting case of T e = T { 



k v {T e = = (1 + Z)k 
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giving that the enhancement factor due to electro-diffusion is (1 + Z). This factor is familiar 
from the well-known ambipolar enhancement for ion diffusion with respect to the laboratory 
frame. However, here the impact is on relative diffusion of two distinct ion species. 

In summary, representing the electric field effect on inter-ion-species diffusion as a modi- 
fication to the conventional baro-diffusion coefficient is only possible when these ion species 
are in the same charge state [Z\ = Zh). Moreover, even in such a case, the electron and ion 
temperatures need to be evolved separately for this effect to be properly accounted for. 
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Appendix A: Evaluating electro-diffusion ratio in the Zel'dovich-Raizer fashion 

In Ref. [5] expression (5) for the baro-diffusion ratio is obtained by utilizing a general 
formula giving k p in terms of the specific volume and chemical potential derivatives over 
concentration. Instead, Zel'dovich and Raizer [6] notice that once k p is known to be a 
thermodynamic quantity, the answer found in some special case should also work for all other 
cases. In particular, evaluating k p can be simplified by considering a globally equilibrated 
system. Indeed, the flux, as well as the temperature gradient, is then equal to zero and the 
baro-diffusion ratio can be obtained by balancing the Vc and V logp terms on the right side 
of Eq. (4). By writing explicit expressions for the densities of the mix components in the 
uniform gravitational field Eq. (5) can then be recovered. Below, we apply this idea to a 
plasma with two sorts of ions. 

First, we recall Eq. (6) and set VTj = = to obtain 




dc 



d log pi eks d& 



0. 



(Al) 



doc ' -L j doc 
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where k p is readily provided by Eq. (5). Next, we notice that for a plasma equilibrated in 
the uniform gravitational field, the light and heavy ion density profiles are given by 

n hh = n m>h0 exp ( - — J , (A2) 

where n^M) are the species number densities at x = and $(x = 0) = 0. Now, the first two 
terms on the left side of Eq. (Al) need to be evaluated with the help of Eq. (A2). 
To do so, we observe that c = m^i /(mini + m^n^ to write 



dc 2 d /1\ 2 d /?Tt/j,n^\ 
dx dx^-c' cfaA mini ' 

Next, we insert Eq. (A2) into the right side of Eq. (A3) to find 



(A3) 



^ = -c(l - c) [-iL(m fc - mi) + ^(Z h ~ Zi)} . (A4) 

Finally, the total ion pressure gradient is calculated along the same lines to obtain 

dlogpi g (mini + m h n h ) eE (Z x n x + Z h n h ) 

— = -T, ,, + ,,, + T t n l + n h ' (A5) 

Then, by inserting Eqs. (5), (A4) and (A5) into Eq. (Al) and solving it for previously 
obtained result (33) is reproduced. 
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